arXiv: 1503.01114v2 [hep-ph] 4Jun2015 


Prepared for submission to JHEP 


CALT-TH-2015-009 


Complementarity of LHC and EDMs for Exploring 
Higgs CP Violation 


Chien-Yi Chen,“ S. Dawson a and Yue Zhang 6 

a Department of Physics, Brookhaven National Laboratory, 

Upton, New York, 11973 

b Walter Burke Institute for Theoretical Physics, California Institute of Technology, 

Pasadena, CA 91125 

E-mail: cychen@bnl.gov, dawson@bnl.gov, yuezhang@caltech.edu 

Abstract: We analyze the constraints on a CP-violating, flavor conserving, two Higgs dou¬ 
blet model from the measurements of Higgs properties and from the search for heavy Higgs 
bosons at LHC, and show that the stronger limits typically come from the heavy Higgs search 
channels. The limits on CP violation arising from the Higgs sector measurements are com¬ 
plementary to those from EDM measurements. Combining all current constraints from low 
energy to colliders, we set generic upper bounds on the CP violating angle which parametrizes 
the CP odd component in the 126 GeV Higgs boson. 
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1 Introduction 

Now that the 126 GeV Higgs boson has been discovered [1, 2], the exploration of its properties 
is the focus of LHC phenomenology. The current measurements of Higgs production and 
decay rates are consistent with the Standard Model (SM) predictions at the ~ 10 — 20% level, 
leaving open the possibility that there is additional physics in the Higgs sector. One attractive 
alternative to the SM is the two Higgs doublet model (2HDM), which has 5 Higgs bosons, 
allowing for new phenomena in the Higgs sector [3]. The couplings of the Higgs bosons to 
fermions and gauge bosons in the CP conserving 2HDM depend on 2 parameters: a, which 
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describes the mixing in the neutral Higgs boson sector, and tan/3, the ratio of Higgs vacuum 
expectation values. Measurements of Higgs coupling properties in the CP conserving limit 
require that the 2HDM be close to the alignment limit, j3 — a ~ ^ [4-7]. 

The SM explains CP violation through the CKM mixing matrix, which is sufficient to 
account for observed CP violation in the B and K systems. However, it is insufficient to 
explain the excess of matter over anti-matter in the universe, suggesting that there may be 
further sources of CP violation [8, 9]. The 2HDM offers the possibility for a new source 
of CP violation beyond the CKM matrix and QCD 6 term. In such a scenario, the 126 
GeV Higgs boson can be a mixture of CP even and CP odd states [10-14] . The LHC 
data has already excluded the case that the 126 Higgs is a pure CP odd scalar [15, 16], but 
the constraints on its CP odd mixture are still rather weak. There have been proposals of 
new techniques to directly measure the Higgs CP mixture in future colliders [17-24]. The 
parameters of the CP violating version of the 2HDM receive complementary limits from 
LHC Higgs coupling measurements and from low energy measurements such as electric dipole 
moments (EDMs). The measurements of Higgs couplings do not put a strong constraint on 
the CP violating phase, especially in the alignment limit [10], and the strongest limits come 
from EDMs [8, 10, 11, 25-27]. 

CP violation in the Higgs sector has been studied extensively in the MSSM limit of the 
2HDM [28, 29]. The MSSM contains many sources of CP violation from the soft SUSY 
breaking terms in the effective Lagrangian [30]. The primary restriction on this type of CP 
violation arises from the requirement that the lightest Higgs boson have a mass near 126 
GeV [31]. Analogous limits to those obtained in this work from Higgs couplings, heavy Higgs 
searches, and EDMs can be found in the MSSM [32, 33]. 

We consider a CP violating 2HDM scenario which has a softly broken Z 2 symmetry 
which avoids large flavor changing neutral currents from Higgs exchange, but allows for new 
CP violation from the scalar potential. We further allow the Higgs couplings to have small 
deviations from the alignment limit. In this work, we consider the additional constraints 
on the parameters of the theory arising from the search for heavy Higgs bosons. In the 
CP conserving 2HDM, the search for heavy Higgs bosons significantly restricts the allowed 
parameter space for small tan /3 [34, 35] and this remains true in the CP violating case. In 
the context of the 2HDMs, if there is significant CP violation, the heavy Higgs boson masses 
cannot be too heavy and in some regions of parameter space the LHC heavy Higgs searches 
can place the leading constraint on CP violation. 

In Section 2, we review the CP violating 2HDM and predictions for Higgs boson produc¬ 
tion and decay within this class of models. Limits from heavy Higgs searches are discussed 
in Section 3 and compared with low energy limits from the electron EDM. We have also 
updated the results of Refs. [10, 36-38] for the limits on the CP violating parameters from 
Higgs coupling fits. Finally, Section 5 contains a concluding discussion of the complementary 
limits on CP violating 2HDMs from Higgs coupling fits, heavy Higgs searches, EDMs, the 
oblique parameters, and g — 2. 
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2 Two Higgs Doublet Models and CP Violation 

In this section we review the 2HDMs considered in this study. 


2.1 Scalar Potential with Two Higgs Doublets 

The most general two Higgs doublet potential which breaks SU(2)l x 17(1) to U(1)em is, 


1 


V(<j)i,<t>2) = [mil ( 0 i 0 i) + (™i 2 ( 0 i 02 ) + h.c.J + 

+y ( 0 i 0 i ) 2 + y (020a ) 2 + A 3 (0i0i)(0 2 02) + A 4 (0i02)(0 2 0i) 


t. 


At, 


1 

+ 2 


A5(01^2) 2 + A 6 (0}0 2 )(0[0i) + A 7 (0j0 2 )(0l,0 2 ) + h.C. 


At, 


At 


t. 


At, 


( 2 . 1 ) 


The potential of Eq. (2.1) leads to tree level flavor changing neutral currents, which can be 
avoided by imposing a Z 2 symmetry under which, 


01 ->■ -01 02 -> 02 ■ 


( 2 . 2 ) 


Eq. (2.2) implies A6 = A 7 = 0, while a non-zero m 42 softly breaks the Z 2 symmetry of 
Eq. (2.2). 

After electroweak symmetry breaking, the Higgs doublets in unitary gauge can be written 
as, 


01 = 


- sin 017+ 

^ (v cos 0 + H® — i sin 0A°) 


= ^ 


02 = e 


cos 


017+ 


^(i;sin 0 + 


H 2 + i cos 0A°) 


(2.3) 


where tan0 = v^fvi , v = \J\ v\ | 2 + | v 2 | 2 = 246 GeV and 17+ is the physical charged Higgs 
with mass m H +. We are free to redefine fields and go to a basis where £ = 0. In general 
there are 2 independent phases and the imaginary parts of mi 2 and A 5 lead to mixing in the 
neutral Higgs sector between and A 0 , and that is the source of CP violation. 

The mixing among the three neutral scalars can be parametrized by an orthogonal matrix 

R, 

^ SaCaib CotCatb ^a b 

7* ^a^a b ^a c C a C ctc S a C OLc C a S ab S ac C Cib S 0ic 

\SaSot b C ac T C a S ac S a S ac C OL S OLb C GLc C OLb C OLc 
where s a = sin a, etc and 

TT TV TV TT . 

~2 <ab ~2 ~ 2 ~ ac ~ 2 (25) 

The physical mass eigenstates are then defined as (hi, /i 2 , h%) T = R(H®, 70?, A°) T . In the CP 
conserving version of the 2HDM, = a c = 0, R is block diagonal, and h\ and /i 2 have no 
pseudoscalar component. 
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2.2 Neutral Scalar Interactions 


For simplicity, we focus on the 2HDMs where the Yukawa sector has a Z 2 symmetry and <p\ 
and (j )2 each only gives mass to up or down type fermions. This is sufficient to suppress tree- 
level flavor changing processes mediated by the neutral Higgs scalars. For the 3rd generation 
(and suppressing CKM mixing), 


C = 


cos amt\yr 

— 77— - 

sin p v 1 

cos amtW 
=—5 — + 

sin p V 


cos a m b \ j- 

-T— o — )Ql ( P20r + h.C. 
sin p V 

sin a m b \ y- 

- a — )Ql ( PiOr + h.c. 

cos p V 


Type I 


Type II, 


( 2 . 6 ) 


where Q ^ = (f^, &l). In both cases, we assume that the charged lepton Yukawa coupling has 
the same form as that of the charge —1/3 quarks. Under the Z 2 symmetry, QlRr ,^ 1 2 are 
always even, <j>i is always odd, and Br is even (odd) in Type I (II) models. 

From this we can derive the couplings between neutral Higgs bosons and the fermions 
and gauge bosons. As a general parametrization we take, 




1=1 


-m f (c f ,Jf + Cf'ifiiBf) ^ + {2a i M^ v W lx W> x + a i M 2 z Z t ^ 


(2.7) 


When CfjCfj / 0 or / 0, the mass eigenstate hi couples to both CP even and CP odd 
operators, so the theory violates CP. The coefficients Cfj, Cf^ and a* can be derived from tan f3 
and the elements of the matrix R defined above. An appealing feature is that all couplings 



Ct,i 

Cb,i — 

Ct,i 

— &T,i 

cii 

Type I 

Rii! sin (3 

Rii! sin/3 

-Rii cot /3 

Ra cot /3 

Ri 2 sin /3 + Rn cos /3 

Type II 

Rii! sin/3 

Rii/ cos /3 

-Ra cot /3 

—Ra tan (3 

Ri 2 sin /3 + Rn cos /3 


Table 1 . Fermion and gauge boson couplings to Higgs mass eigenstates. 


in Table 1 depend on only four parameters, a, a c and tan/3. It is worth noting that the 
couplings of the light Higgs boson hi to the gauge bosons and fermions are independent of 
a c - Fits to the CP conserving 2HDM suggest that the couplings are close to the alignment 
limit, j3 — a ~ implying that h\ has couplings very close to the SM predictions. In our 
numerical studies, we will allow small deviations from the alignment limit. 

2.3 CP Violation Implies a Non-Decoupled Heavy Higgs Sector 

In general, the imposed Z 2 symmetry in the Yukawa sector is not preserved by renormaliza¬ 
tion. The hard breaking Ag, A 7 terms from the Higgs potential will induce couplings of 4> 1 , 4 >2 
to both up and down type quarks. This does not reintroduce any tree level flavor changing 
effects because the induced Yukawa matrices are still aligned with the corresponding fermion 
mass matrices. Motivated by this, a convenient choice is to forbid the Ag, A 7 terms. In this 
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case, the model has an approximate Z 2 symmetry, which is only softly broken by the m\ 2 
term. 

For the approximate Z 2 symmetric model, all of the potential parameters can be solved 
for in terms of the following parameters: 

• The scalar masses, m^, 'm/, 2 , nrih 3 and m H ± 

• The neutral scalar mixing angles, a, otb, a c 


• The ratio of vev’s, tan /3 

• One potential parameter, Re(mf 2 ), or v = Re(mi 2 ) 2 /(u 2 sin 2/3) , 


giving 9 physical parameters. The v parameter controls the decoupling limit, he., when 
Re(mi 2) 2 approaches infinity, the masses of h. 2 , h% and H± also go to infinity. 

The explicit solution for the parameters of the scalar potential was found in Ref. [10], and 
is listed below in Appendix A. The imaginary part of A 5 , which is a source of CP violation, 
is given by, 


2 cos afe / 2 2 \ 

ImAs = —^-- mi, - miJcosasma c cosa c 

v z sm pi J 


+(jn? hl — m\ 2 sin 2 a c — m 2 3 cos 2 a c ) sin a sin 


( 2 . 8 ) 


An important point here is that, in order for the 126 GeV Higgs boson to have CP violating 
couplings, the heavy Higgs states must not decouple. Otherwise, the two Higgs doublet model 
will return to the SM limit. This is actually our main motivation for studying the bounds on 
the non-decoupled heavy Higgs. 

Clearly, when the scalars /i 2,3 are much heavier than the electroweak scale, and mh 2 — 
m h3 = m H + > m h[ , 


sin 2 q '5 |~ 


ImA 5 | v 2 


m 


H+ 


sin /3 
sin a 


(2.9) 


The unitarity bound on ImAs, ImAs < 4tt, sets the largest allowed CP violating mixing angle 
Ub- This implies that for an 0(1) sin ct;, to be theoretically accessible, the heavy scalars h 2, /13 
and H± must be not far above the electroweak scale. In general, for nonzero ctb, the masses 
of the heavy scalars should satisfy 


m H+ Ip 870 GeV x -y/| ImAs | /(4vr)^/| sin / 0/(sinasin2a6) | . (2-10) 


A similar conclusion holds when one goes beyond the approximate Z 2 symmetry by including 
the Ag, A 7 terms. 
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2.4 Beyond Approximate Z 2 Symmetry 


For the approximate Z 2 symmetric model, there is a further theoretical constraint on the 
physical parameters resulting from the minimization of the potential. This constraint is given 
in Eq. (A. 10) and can be transformed into a quadratic equation for tana c . The condition for 
a c to have a real solution is 


sin 2 ajj < 


(m 2 h3 -m 2 h2 ) 2 cot 2 (a + f3) 
4(m 2 h2 - m 2 hi )(m 2 h3 - m 2 hi ) 


= sin 2 


_ max 

a b 


( 2 . 11 ) 


When Eq. (2.11) is satisfied, the solutions for a c are, 


OLr = 


a c , a + f3 < 0 
a+, a + (3 > 0 


tan a± = 


=F| sin a™ 


sin 2 aP ax — sin 2 a b 


sm ab 


I m 


h 3 


— m 


hi 


m 


hi 


— m 


hi 


( 2 . 12 ) 


Eq. (2.11) implies an additional theoretical upper bound on the CP violating angle ab, 
when the other parameters are fixed. In practice, we sometimes find this bound can be 
stronger than all the experimental limits. However, this is only a bound because of theoret¬ 
ical prejudice. In fact, it can be removed with a minimal step beyond the approximate Z 2 
symmetric case by introducing a A 7 term, with A 7 being purely imaginary. In this case, the 
bound Eq. (2.11) no longer exists, a c becomes a free parameter, and ImA 7 can in turn be 
solved for as, 


ImA7 = 


2 cos ab 
v 2 tan 2 /3 


(m 2 h , — rn'b 2 ) sin a c cos a , 


cos(a + /3) 
cos 2 f3 


,/2-2 , 2 2 2 \ • 

+(m . h2 sm a c + m h3 cos a c — m hl ) sm ab 


sin(a + /3) 
cos 2 (3 


(2.13) 


Although introducing hard Z 2 breaking (A 6,7 7 ^ 0) makes the Yukawa structure in Eq. (2.6) 
unnatural, one might argue it is accidentally the case at the electroweak scale. 1 In the 
phenomenological study in the next section, we will give the results for both the approximate 
Z 2 case, and the minimal extension as discussed in this subsection. 


2.5 Production and Decay of the Heavy Higgs at LHC 
2.5.1 Production 

The dominant heavy Higgs boson production channels relevant to this study are gluon fusion, 
gg -» /i 2 , 3 , vector boson fusion, qq —> qqh 2 , 3 , and production in association with bottom 
quarks, gg/qq h- 2 , 366 . In the 2HDM we explore, the interactions between the heavy neutral 

1 We are aware that allowing Zi breaking terms in the Yukawa sector can introduce additional sources of 

CP violation. The price for this is introducing tree level flavor changing effects at the same time, and some 
flavor alignment mechanism must be resorted to [39-41]. We do not consider such a possibility, but focus on 
CP violation only from the Higgs sector in this work. 
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Higgs bosons and the SM fermions and the W, Z gauge bosons are simply rescaled from those 
of a SM-like Higgs boson, H$m, by a factor given in Table 1. Therefore, it is convenient to 
take the SM-like Higgs cross sections, and rescale them with these factors and the appropriate 
form factors. The LHC production cross sections for a heavy SM-like Higgs boson have been 
calculated by the LHC Higgs Cross Section Working Group and given in [42, 43]. 

For the gluon fusion process, we calculate the ratio of the heavy Higgs boson production 
cross section in a CP violating 2HDM to that of a SM-like Higgs with the same mass. At 
one-loop, 


EL. = 


p(gg -t L) 

ct^ /2 (rl) + c feii Af / 2 (r b l ) 

2 

+ 

Ct,iAf /2 (Ti) + c b}i Af /2 (T l b ) 

°(gg -t Hsm) 

A i /2 ( T t) + A l/2^ T b) 

2 


,(2.14) 


where r) = m)],/(4m 2 ) and i = 1,2, 3, / = t, b. The form factors A^ 2 , Ay 2 are given by 



fir) 


2(r + (r- 1 )/(t))t 2 , 
2 /( r ) r_1 , 


arcsirr 


(vV), 



( l+\/l — T 1 \ 

VW1-T-V 



T < 1 
T > 1 


(2.15) 

(2.16) 

(2.17) 


For vector boson fusion, the ratio of the heavy Higgs production cross section in a 2HDM 
to that of a SM-like Higgs with the same mass is simply 


R'vrf — 


a(qq -)• qqhj) 
a(qq qqH S M) 


(a,) 2 . 


(2.18) 


For /i 2 i 3&6 associated production, we take the NLO cross section for SM-like Higgs boson 
production in the 4 flavor number scheme from Ref. [44, 45]. There the cross section contains 
two pieces, one is proportional to y 2 , and the other proportional to y b yt from interference. 
We rescale these results with the heavy Higgs-fermion couplings in a 2HDM, 


a(bb ->• hi) = (c bti ) 2 a b (m hi ) + (m hi ) + (c b ,i) 2 ° b {m hi ) + ct : iC bti af(m hi ) , (2.19) 

where is the cross section for gg —> bbhi where the Higgs couples to the b quarks, af 1 
is the interference between diagrams contributing to gg — > bbhi where the Higgs couples to 
the b and the t quark, and af are the corresponding contributions from the pseudoscalar 
couplings to the b and t quarks given in Eq. (2.7). Results in the 5 flavor number scheme [46] 
are quite similar and do not affect our conclusions. 


2.5.2 Decays 

The heavy neutral scalar to electroweak gauge boson decay rates are 


r (hi ->■ vv) = (oi) : 


Gfiti' 


hi 


lev^vr 


VI- 


4M 2 \ 

m l) 


1/2 



( 2 . 20 ) 
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where V = W, Z and d\y = 2, 6 Z = 1, and i = 2,3. We note that in the alignment limit, 
r(h 2 ,3 —> VV) = 0 when sina^ = sina c = 0. These channels open up with non-zero CP 
violation. The decay rates to SM fermions are 


r (hi ->• If) 


0 Of,i ) 2 + ( Cf,i ) 2 


N c G F m 2 m h . 
Ay/2i r 



( 2 . 21 ) 


where N c = 3 for quarks and 1 for charged leptons. 

The heavy scalars can also decay to a pair of gluons via a loop of top or bottom quarks, 
and the rates are 


r (hi ->• gg) = 


a sGFmj. 

64\/2vr 3 


3 r 


c t,i^i/2i T t) + c b,i^l/2i T b) + + 4^f/2( T fe) 




.( 2 . 22 ) 


Clearly a decay rate is a CP even quantity. Thus, in all the above decay rates, the CP even 
coefficient c l j and the CP odd one always contribute incoherently. 

In our study, we are also interested in the heavy neutral scalars, h 2 , h. 3 , decaying into the 
Z boson and the 126 GeV Higgs boson, 


r (hi Zhi) = J (™ 2 h . - (m hl + M z ) 2 ) (m 2 h . - (m hl - M z ) 


-( 2 ™ 2 h . + 2 m 2 hl - Mf) + (m£. - m 2 hl ) 2 

z 


(2.23) 


where g lz \ = (e/ sin29w) [(- sin f3R u + cos/3i?i 2 )i?j 3 - (—sin/3Rn + cos/3i?,; 2 )i?i 3 ]. 

We have also calculated the decay rate of hi —>• 2hi from the Higgs self-interactions. The 
decay rate is 


Y{hi ->• hihi ) = 


2irm hi 


' 1 - 


4m? 


m 


fei 

l 

hi 


(2.24) 


where < 7 * 11 , (?' = 2,3) are defined in Appendix B. 

To get the branching ratios, we calculate the total width of the heavy Higgs 2 , 


r totihi) = T{hi ->■ w+w~) + T{hi -> zz) + r(h, ->■ tt) + r(h* -> hi) 

+ r(/ij —> t+t ) + r(hj —> gg) + r(hj —> Zh\) + r(hj —> hi hi) . (2.25) 

Finally, for each channel, the ratio of signal strengths in the 2HDM to the counterpart in the 
SM is given by, 


xx 


hi 


(a\C 7 + 8 ) x Br(h, : -> XX) 

(crf M £ 7 + crf M £ 8 ) x Br SM (hi ->• XX) ’ 


The rate hi —> 77 for i = 2,3 is always small and can be neglected here. 


(2.26) 
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where, for example, the production cross sections are given by 

°7 = a 99-7R l gg + ^VBFjRyBF + a VH,7^VH > (2.27) 

(Jggj is the gluon fusion cross section from Ref. [42, 43] for a SM Higgs boson with a mass of 
to/,,., and £ 7 ^ are the luminosities used in the experimental analysis. With this quantity, we 
are able to reinterpret the constraints on a heavy SM-like Higgs boson for the heavy neutral 
scalars in the 2HDM. 


2.6 CP violation and Heavy Higgs Signal Rates 


At this point, it is useful to gain some intuition about the impact of CP violation on the 
heavy Higgs to gauge boson decay channels, hi —> VV and /q —> Zh\ with (i = 2,3). It is 
convenient to redefine the Higgs doublets and go to a basis where only one doublet, called 
4>'ii gets the 246 GeV vev, while the other (j )' 2 has no vev [13, 67]. 

We start from a special point in the parameter space where the lightest Higgs, hi, has 
exactly the same couplings does the SM Higgs boson. This corresponds to having the mixing 
angles in Eq. (2.4) satisfy otb = a c = 0, and f3 — a = n/2. The Higgs sector preserves CP 
invariance at this point. In this case, h\ is the excitation arising from <f>\ defined above, while 
Tj- 2,3 are excitations from </> 2 . As a result, the decay rates hi -7 VV and hi —> Zh± both vanish 
for i = 2,3. It is worth noticing that this special point can be approached without going to 
the real decoupling limit by sending the second doublet mass to infinity. 

Next, we turn on CP violation by making otb = 0.5, but still keep a c = 0. Here we discuss 
an example by fixing tan f3 = 20 (in the basis of {0i , 02 } given in Eq. (2.3)) and vary the 
angle a, or the quantity cos(/3 — a). We also choose the heavy neutral scalar masses to be 
mh 2 = 400 GeV and mg l3 = 450 GeV. In Fig. 1, we plot the gluon fusion production cross 
section and the gauge boson branching ratios of h. 2,3 as a function of cos (/3 — a). There are 
several suppressed regions which can be understood from Table 1. In the case a c = 0, we 
have in the Type-I model, 


sin a _ 

Ct,2 = c& 2 =-;— Ct ,2 = — 0, 2 = 0 , 

sin p 

COS OL 

ct ,3 = c-b 3 = —;—— sin ab, ct 3 = —q, 3 = — cos otb cot /3 . 
snip 


(2.28) 


First, the gluon fusion production cross section for /12 via a top or bottom loop vanishes at 
a = 0 (near cos(/3 — a) ~ 0). In the example we describe here, f3 = arctan(20) is close to 
7t/2, and ctp = vanishes at a = 0. Second, at a = ±7t/2, (near cos(/3 — a) ~ ±1), the 
couplings ct ,3 and Cb ,3 vanish. As a result, the production cross section for /13 is suppressed 
because q .3 and are both suppressed by cot/3 = 1/20 in this case. On the other hand, 
the gauge boson decays of h 2,3 are directly controlled by (3 — a. We list the relevant couplings 
here, again for a c = 0 , 


02 = — cos (/? 
a 3 = - sin(/3 


a): 92zi — — 


sin 2 9w 


a)sinafe, g^zi = ~ 


sin 29 w 


sin(/3 — a) sin op, 
cos (P — a) , 


(2.29) 
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1 


hj h 3 



Figure 1. An example showing the impact of a non-zero CP violating angle, ab = 0.5, and the 
deviation from alignment (parameterized by cos(/3 — a)) on the heavy Higgs production from gluon 
fusion at y/s = 8 TeV (left panels) and their decays (right panel) in h n — ► VV (red, solid) and 
hi —► Zh\ (blue, dashed) channels. We have fixed the other parameters to be tan/3 = 20, a c = 0, 
m-h 2 = 400 GeV, mh 3 = 450 GeV and u = 1. 


where gi z \ (i = 2, 3) is the coupling between hi-Z-h\ defined below Eq. (2.23). These make it 
manifest why the heavy Higgs to gauge boson decay channels are sensitive both to a deviation 
from the alignment limit and to CP violation. Clearly, when cos(/3 — a) = ±1, the decay 
rates /13 — > VV and /12 -A Zh\ vanish, while when cos (j3 — a) = 0, the decay rates —> VV 
and /13 —> Zh\ vanish. For the case of /12 decay, the branching ratios are more suppressed 
because the decay h 2 —> h\h\ dominates in most of the parameter space. Therefore, the most 
important constraints come from the /13 — > VV and /13 — > Zh± channels. 

Combining Eqs. (2.28) and (2.29), we find the /13 — > VV signal rate (production cross 
section x decay branching ratio) is peaked at cos(/3 — a) = 0, while /13 —> Zh\ vanishes at 
both cos (/3 — a) = 0, ±1, and is peaked in between. With these facts, one can understand the 
yellow and orange regions in the upper right panel of Fig. 11. One can also follow a similar 
analysis in order to understand the generic features in the other plots. 
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3 Results 


In this section, we describe our method to obtain constraints from heavy Higgs searches at 
the LHC, and show the numerical results in a series of figures. 

In the presence of CP violation, all of the three neutral scalars mix together, and we 
fix the lightest scalar, h\, to be the 126 GeV scalar already discovered at the LHC. As 
discussed in the previous sections, the heavy Higgs to gauge boson decay channels, including 
^2,3 H ► WW/ZZ and h 2,3 — > Zh\ — > l + l~bb, are not only sensitive to deviations from the 
alignment limit (/3 — a = 7 r/ 2 ), but also to the presence of CP violation ( ab,a c 7 ^ 0). We use 
the production and decay rates calculated in Sec. 2.5 to obtain the 2HDM predictions for the 
heavy Higgs signal strength in these two channels. Then we compare these theory predictions 
to the results from the 7 and 8 TeV running of the LHC. For the heavy Higgs search data, 
we use limits for masses up to a TeV from the h 2,3 —> WW/ZZ channel [48, 49] and from the 
h- 2,3 —> Zh\ — > l + l~bb(r + T~ ) channel [50, 51]. 

We also take into account the /i 2 ; 3 —> t + t~ channel [52], which gives constraints for heavy 
Higgs masses up to a TeV and is relevant in the Type-II model in the large tan/3 case [53]. 
The experimental results are given as correlated fits to a(bb — > h)BR(h —> t + t~ ) versus 
cr(gg —> h)BR(h —> t + t~ ), which are shown in Fig. 8 of Ref. [52]. We interpret these bounds 
as bounds on the production of heavy Higgs states in the CP violating 2HDM. For the case 
nrih 2 = rrih 3 = 300 GeV, the limits are shown in Fig. 2 3 . In addition to a limit for tan/3 > 30, 
there is another limit around tan/3 ~ 1. This lower limit arises because the masses of h. 2,3 
are below the tt threshold and hence the dominant decay channel is through h 2,3 —> bb and 
fr- 2,3 —» t + t~. For mh 2 = m/i 3 = 500 GeV the bound from /i 2,3 —> t + t~ gives no constraint 
for tan/3 < 50. In general, we find that the bounds from h. 2,3 —> t + t~ are always weaker 
than those from the coupling measurement of the light Higgs or EDMs. As a result, we will 
not include them in the following plots. 

The most up-to-date 126 GeV Higgs coupling data are given in Table 2, normalized to 
the appropriate luminosities. They are used to constrain the theoretical predictions for the 
signal rates of h\, from Sec. 2.5. We take the SM cross sections from the LHC Higgs Cross 
Section Working Group [43]. We have performed a x 2 analysis using the results listed in 
Table 2. 

In Figs. 3 to 12, we show the limits derived from heavy Higgs searches and the light (126 
GeV) Higgs data, together with those from the low energy electron and neutron EDMs. For 
the EDM constraints, we use the results of Ref. [10]. 

In these numerical results, we fix the heavy Higgs masses and the parameter v = 1. The 
CP violating angle a c is fixed in the approximate Z 2 symmetric model by Eq. (2.12). On 
the other hand, for the extended model without an approximate Z 2 symmetry, a c is a free 
parameter. We also note that varying the parameter v between 0 and 1 only leads to slight 
changes to our results. The constraints are shown in the sin ab versus tan /3 plane, while 
varying a and a c . We consider both the alignment limit with a = /3 — 7 t/2 and cases when 

li We assume that there is no interference between the /12 and hi resonances. 
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Type-II 2HDM, alignment limit, m hl 3 =300 GeV 



Figure 2. Heavy Higgs search constraints on the Type-II 2HDM with approximate Zi symmetry and 
rrih 2 = rrih 3 = 300 GeV. The colored regions are excluded by the search for h 2 ,3 —> t + t~ . 


Table 2. ATLAS and CMS Higgs Coupling Measurements. 


Channel 

VCMS 

Ref. 

HATLAS 

Ref. 

Hww 

0.83 ±0.21 

[54] 

1 09 +0 ' 23 
i ' uy - 0.21 

[55] 

Vzz 

1.0 ± 0.29 

[54] 

1 44 + 0 - 40 

[56] 

M 77 

1.13 ±0.24 

[54] 

1.17 ±0.27 

[57] 

l^bb 

0.93 ±0.49 

[54] 

0.5 ±0.4 

[58] 

Htt 

0.91 ±0.27 

[54] 

1.4 ±0.4 

[59] 


there are small deviations from alignment, cos(/3 — a) = ±A. The 126 GeV Higgs data put 
upper bounds on A for fixed values of tan/3. For the Type -I model, we consider A = 0.1, 
while for the Type- II model, the light Higgs coupling data constraint is stronger at large 
tan/3, so we take A = 0.02 [4]. ATLAS has also limited the parameters of the 2HDM by 
directly searching for the heavier neutral Higgs boson, but these limits are not competitive 
with the Higgs coupling data for the heavy h 2,3 masses that we consider [47] . 

3.1 Limits from Heavy Higgs Searches in Approximate Z 2 Symmetric Models 

We first discuss the models with an approximate Z 2 symmetry. Fig. 3 shows the limits on the 
CP violating parameter, a &, as a function of tan /3 in the Type-I model. In each panel, the gray 
area marked “theory inaccessible” has no real solution for a c from Ecp (2.12). The left panel 
assumes the alignment limit, /3 — a ~ 7 r/ 2 , while the right panel allows for a small deviation 
from the alignment limit 4 . The orange area is excluded by the heavy Higgs search channel 
Tj- 2,3 -a Zh\ —> l + l~bb , while the yellow area is excluded by the channel /i 2,3 —> WW/ZZ. It is 
clear that the limits become quite stringent away from the alignment limit. For comparison, 

4 The results are similar for negative A. 
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Type-I w. approximate Z 2 : m/, 2 =400GeV, =450GeV, m H ± =420GeV, v=l 

cos(/J-ff)=0.1 
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Figure 3. Heavy Higgs search constraints on the Type-I 2HDM with approximate Z 2 symmetry, using 
the /i2,3 —> WW/ZZ (yellow) and ^2,3 —> Zh\ —> l + l~bb (orange) channels. These constraints are 
presented in the sin at versus tan j3 parameter space and colored regions are excluded. The left panel 
is for the alignment limit with a = /3 — 7t/2, while the right panel shows the case with a deviation 
from that limit. Also shown in blue are the electron EDM excluded regions. In these plots, we have 
chosen the heavy scalar masses to be rrih 2 = 400 GeV, rrih 3 = 450 GeV, m H + = 420 GeV, and the 
model parameter v = 1. The other mixing angle a c is a dependent quantity fixed by Eq. (2.12). In 
the gray region, there is no real solution for a c . 

Type-I w. approximate Z 2 : m/, 2 =550GeV, =600GeV, m H ± =620GeV, v=l 

f3-a=n/2 cos(J3-a)=0A 




Figure 4. Similar to Fig. 3, but with heavy scalar masses m/, 2 = 550 GeV, rrih 3 = 600 GeV, m H + = 
620 GeV. In the right panel, the red region is excluded by the 126 GeV Higgs data applied to /q. 

we include the results of Ref. [10] for the limits from the electron EDM (eEDM, the blue 
shaded regions are excluded). In all cases, the EDM limit and the heavy Higgs searches 
exclude complementary regions. The masses of the heavy Higgs are increased to around 600 
GeV in Fig. 4. In this case, the limits from heavy Higgs searches become much weaker, with 


13 - 






















Type-II w. approximate Z 2 : m/, 2 =400GeV, m/i 3 =450GeV, m H ± =420GeV, v= 1 
fi-a=nl 2 cos(/3-ar)=0.02 




Figure 5. Heavy Higgs search constraints on the Type-II 2HDM with approximate Z 2 symmetry. 
The Higgs sector parameters are chosen to be the same as those in Fig. 3. The colored regions are 
excluded by searches for /i 2 ,3 —> WW/ZZ (yellow), h 2 ,3 —► Zh\ —> l + l~bb (orange) channels, 126 GeV 
Higgs coupling data (red), electron EDM measurements(blue), and neutron EDM limits(green). The 
gray region is again theoretically excluded because it contains no real solution for a c . 


Type-II w. approximate Z 2 : m/, 2 =550GeV, m/, 3 =600GcV, m H± =620GcV, v= 1 
p-a=n/2 cos()8-q')=0.02 




Figure 6. Similar to Fig. 5, but with heavy scalar masses to^ 2 = 550 GeV, m ^ 3 = 600 GeV, m H ± = 
620 GeV. 

the dominant excluded region coming from the eEDM searches. The mass splitting between 
the heavy masses is restricted by limits on the oblique parameters, which is discussed in Sec. 
4.1. 

Fig. 5 shows the limits on versus tan j3 in the Type-II model. Away from the alignment 
limit (the right panel), there is a significant exclusion region for tan (5 > 10 from the 126 GeV 
Higgs parameter measurements. Around tan/3 ~ 1, the electron EDM constraint vanishes 
due to a cancellation among the Barr-Zee diagrams as pointed out in Ref. [ 8 ]. We find the 
heavy Higgs searches from the gauge boson decay channels h. 2,3 —> Zh\ and /i 2,3 —> WW/ZZ 
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are extremely useful and close the window of large values of sina;, ~ 0(1) in all cases. As 
the mass of the heavy particles is increased in Fig. 6, the region excluded by the heavy Higgs 
searches shrinks, with again the dominant exclusion coming from the eEDM and neutron 
EDM (nEDM, the green regions are excluded). It is worth pointing out that the neutron 
EDM excluded regions are shown using the central values given in [60] , which however involves 
large uncertainties in the evaluation of hadronic matrix elements. In contrast, the heavy Higgs 
searches provide a robust upper limit on the CP violating angle ab- 

3.2 Limits from Heavy Higgs Searches in the Models with no Z 2 Symmetry 

As discussed in Sec. 2.4, if the assumption of an approximate Z 2 symmetry is relaxed, the 
theoretical relationship between a\, and a c can be removed. In this case a c becomes a free 
parameter. This helps to remove the theoretically inaccessible region in Figs. 3-6, and one 
can get a complete view of various constraints in the whole parameter space. 

Figs. 7 and 8 show the constraints in the Type-I model with a c chosen equal to 0 or ab, 
and with two sets of heavy Higgs masses. It is apparent that the dependence on a c is rather 
weak. The results in the Type-II model are shown in Figs. 9 and 10. In Figs. 11 and 12, the 
heavy Higgs search constraints are also displayed in the ab and cos ((5 — a) plane. 

It is also worth re-emphasizing that at low tan/3 ~ 0(1) the 126 GeV Higgs data puts a 
very weak constraint on the CP violating angle a&. This can also be understood from Table 1, 
where the lightest (126 GeV) Higgs couplings to other SM particles near the alignment limit 
are 


aq ~ cos ab , ct,i — (1 + A cot /?) cos ab , 

(1 + A cot /3) cos ab, Type — I 
(1 — A tan /3) cos a;,, Type — II 

where A = cos(/3 — a) and we have kept terms up to first power in A. Clearly for small 
A and tan/3 ~ 1, all CP even couplings are approximately cos ab and all CP odd couplings 
~ ± sin ab- They approach the values in the SM limit when ab —> 0. In the presence of CP 
violation, the gluon fusion production cross section of the light Higgs gets rescaled from the 
SM value by a factor [ 8 , 10], a(gg —> h\)/(J sm(99 —> H$m) — 1 + 1.42 sin 2 afi. The vector 
boson fusion and associated production rates get suppressed by a(VV —> h\) / as\[(VV — > 
Hsm) = cos 2 cq,. The light Higgs to fermion (h\ —> bb,r + r~) decay rates are not affected 
because the CP even and CP odd couplings contribute incoherently, T(h\ —> f f)/TsM.(h -» 
//) = cos 2 ab + sin 2 ab = 1. The Higgs to gauge boson (h\ —> WW* , ZZ*) decay rates get 
suppressed, r(/ii -> W*)/Tsm(/t —> VV*) = cos 2 ab- The light Higgs to diphoton decay 
rate in the presence of CP violation has been given in refs. [ 8 , 10], which in this case can be 
simplified to r(/ii —> 77 )/Tsm (h —> 77 ) ~ 1 — 0.81 sin 2 As a result, the final y 2 of the fit 
for the 126 GeV Higgs data depends on cos 2 ab, and for the SM case y| M = y 2 (cos 2 ab —> 1). 
Because the cos 2 a;;, function is very flat near ab = 0, one can maintain a fit as good as in the 
SM for sizable ab- 



C t A 


Chi - 


— cot /3 sin ab , 

cot /3 sincq,, 

— tan (3 sin ab, 


Type - I 
Type - II 


(3.1) 
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In contrast, the heavy Higgs decay to gauge boson channels (h' 2,3 —> VV and Zh\) are 
more sensitive to a non-zero CP violating angle ab and can place a stronger constraint on it. 
This feature has been discussed in Sec. 2.6. Furthermore, from the figures we notice that at 
low tan /3, the heavy Higgs search constraint is stronger than at large tan /3. This is because 
hi —>• tt, (i = 2,3) is the dominant decay mode and the branching ratio for the gauge boson 
decay modes of hi can be written as 

Br^-^yy or zh (l° w tan/3) ~ Th i~/ vv or Zh _ (3.2) 

hi—ytt 

Eqs. (2.28) and (2.29) tell us that these two rates around the alignment limit are both in- 
sensitiveto variations of tan/3. However, as tan/3 grows to larger than 0(2), the other decay 
channels such as hi —> h\h 2 and hi —> bb larger than hi —> tt, and they are not yet con¬ 
strained by the LHC data. As a result, the gauge boson decay rates of heavy Higgs bosons 
are suppressed in this region. 

Fig. 13 depicts 95% CL constraints in the tan/3 versus cos(/3 — a) plane from heavy Higgs 
searches (black) and from 126 GeV Higgs data (yellow) on the Type-I (first row) and Type- 
II (second row) 2HDMs without approximate Z 2 symmetry. Different curves correspond to 
ctb = 0 (dotted), 0.1 (solid) and 0.5 (dashed), and the other mixing angle a c = 0. For the 
CP conserving case (otb = 0), we found that the bounds are very similar to those studied 
in Refs. [5, 34, 35, 61]. In both Type-I and Type-II models, both heavy and light Higgs 
searches favor regions around the alignment limit cos(/3 — a) = 0. In the Type-II model when 
CP violation is small (bottom left panel), there is another allowed branch corresponding to 
cos(/3+a) ~ 0 [62], but we find the heavy and light Higgs favored regions are inconsistent with 
each other for very large deviations from the alignment limit. In the Type-I model (first row), 
the light Higgs bound only depends on cos(/3 — a) , but is independent of tan/3 in the large 
tan/3 limit. The reason is that in this case the h\ couplings can be approximated as ct,i = 
CbA —t sin(/3 — a) coscq, + 0(1/ tan/3), ct.\ = —c&q = 0(1/ tan/3) and aq = sin(/3 — a) cos ab, 
so their dependence on tan /3 is suppressed. On the other hand, for the Type-II model, the 
couplings C 5 . 1 ,c Ti i and q,. 1 ,<%, 1 are enhanced at large tan/3. This explains why in the Type-II 
model (second row), light Higgs data are more restrictive on the parameter space with large 
tan/3. As a result, for ab = 0.5, the light Higgs data only favors a region with tan/3 < 2 
(see the bottom right panel of Fig. 13). In contrast, we have learned that the heavy Higgs 
search data are more sensitive at small tan (3 and for ab = 0.5 they only allow the region 
where tan (3 > 3, thus there is no region in the parameter space that can be made consistent 
with both light and heavy Higgs results from LHC. Fig. 14 gives results similar to those 
in Fig. 13 but with a different set of mass parameters, rrih 2 = 550 GeV, mh 3 = 600 GeV, 
m H + = 620 GeV. The parameter space becomes less constrained by the heavy Higgs searches 
because the production cross sections are smaller compared to those in Fig. 13. 

From the above results, we can conclude that if the heavy Higgs masses lie below around 
600 GeV, the CP violating phase ab is constrained to be less than around 30% throughout 
the most general parameter space. The regions which allow ab close to this upper bound are 
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tan f3 ~ 1 in the Type-II model, and tan f3 > 20 in the Type-I model without an approximate 
Z 2 symmetry. We have also estimated the future sensitivity of the heavy Higgs search at the 
14 TeV LHC by rescaling the current limits by the square root of expected number of events 
(<7 x C). With 300 (3000) fb -1 data, if the heavy Higgs masses are below 600 GeV and we 
still do not find them, the CP violating angle at, will be constrained to be less than around 
10 %. 

Recall that the angle ab parametrizes the size of CP odd mixture in the 126 GeV Higgs 
boson. The main point of this work is to show that the heavy Higgs search is relevant and 
plays a complimentary role to the other indirect searches, and sometimes it stands at the 
frontier of probing the Higgs boson CP mixture. 
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Type-I w/o Z 2 : =400GeV, m/, 3 =450GeV, m H ± =420GeV, v= 1 




0.5 1 2 5 10 20 50 


0.5 1 2 5 10 20 50 



Figure 7. Heavy Higgs search constraints on the Type-I 2HDM without approximate Z 2 symmetry, 
i.e ., in this case a c is a free parameter which is allowed to vary. The color scheme for the exclusion 
regions is the same as in Figs. 3-6. The first two rows use the same parameters as Fig. 3, and the last 
two rows use the same as Fig. 4. 


Type-I w/o Z 2 : m/, 2 =550GeV, =600GeV, =620GeV, v= 1 




tan ft 


cos(/?—ar)=0.1, ar c =0 
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Figure 8. Similar to Fig. 7, but with heavy Higgs masses rrih 2 = 550 GeV, mu 3 = 600 GeV, m H + = 
620 GeV. 
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Type-11 w/o Z%: m hl =400GeV, m/, 3 =450GeV,m w± =420GeV, v=l 
p-ct=n[2, a c =0 cos(^-a)=0.02, a c =0 






Figure 9. Similar to Fig. 7, but for the Type-II 2HDM without approximate Z 2 symmetry. 


Type-II w/o Zfe: m hl =550GeV, =600GeV, m H± =620GeV, v=l 

j3-a=n/2, a c =0 cos(J3-a)=0.02, a c =0 






Figure 10. Similar to Fig. 8, but for the Type-II 2HDM without approximate Z 2 symmetry. 
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Figure 11. Heavy Higgs search constraints on the Type-I (first row) and Type-II (second row) 
2HDM without approximate Z 2 symmetry, using the /i 2,3 —> WW/ZZ (yellow) and /i 2,3 —► Zh\ —> 
l + l~bb (orange) channels. The heavy scalar masses are fixed to be m^, 2 = 400 GeV, m^ 3 = 450 GeV, 
to#+ = 420 GeV, and the model parameter v = 1. The other mixing angle a c = 0. 
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Figure 12. Similar to Fig. 11, but with heavy scalar masses rrih 2 = 550 GeV, mh 3 = 600 GeV, 
m H + = 620 GeV. 
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Figure 13. Heavy Higgs search (black) and 126 GeV Higgs data (yellow) constraints at 95% CL on 
the Type-I (first row) and Type-II (second row) 2HDM without approximate Z-i symmetry. Different 
curves correspond to otb = 0 (dotted), 0.1 (solid) and 0.5 (dashed). The heavy Higgs curves include the 
combination of constraints from /i2,3 WW/ZZ, h 2,3 —> Zh\ —» l + l~bb and h 2,3 —> t + t~ channels. 
The heavy scalar masses are fixed to be /rp , 2 = 400 GeV, m^ 3 = 450 GeV, m H ± = 420 GeV, and the 
model parameter v = 1. The other mixing angle a c = 0. 
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Figure 14. Similar to Fig. 13, but for heavy Higgs masses rrih 2 = 550 GeV, rrih 3 = 600 GeV, m H ± = 
620 GeV. 
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4 Limits from B Decays, Oblique Parameters, and (g — 2) M 


The CP violating 2HDM is also limited by measurements in B decays, the oblique parameters, 


and (g — 2) M . In Type-II models the charged Higgs mass is restricted by B data to be 
greater than m H + ~ 340GeV for all values of tan/3. In both Type-1 and Type-2 models, 
measurements in the B system prefer tan/3 > 1 [5, 63, 64]. 

4.1 Limits from Electroweak Oblique Parameters 

The allowed parameters are restricted by measurements of the oblique parameters. The 
general results for S,T and U in a 2HDM are given in Refs. [3, 65-67]. In the alignment 
limit, cos ol = sin/3 and sin a = — cos/3, the results simplify considerably, 


jsin 2 a>bF(m 2 H+1 m 2 hl ) + (1 - sin 2 a^sin 2 a c )F(m 2 H+ ,m 2 h2 ) 

+(1 — sin 2 a b cos 2 a c )F{rn 2 H+ ,m 2 h3 ) — cos 2 a c sin 2 ct b F(m % ll , m 2 2 ) 
— sin 2 a c sin 2 a b F(m 2 1 , ) - cos 2 a b F(m \ 2 , m| 3 ) 

+3cos 2 a b F(M 2 z ,m 2 hl ) - F(M^,m 2 hl ) 

+3 sin 2 a c sin 2 aib F(Ml,m 2 h2 ) - F(M^,m 2 h2 ) 

+3 cos 2 a c sin 2 a b [F{M z , m| 3 ) - F(M 2 V , m/j 3 ) 




(4.1) 


where the last line is the subtraction of the SM Higgs contribution evaluated at the reference 
scale, Mjj ref , at which the fit to the data is performed. The function F(x,y) is, 


F(x,y) 

F(x, x ) = 0 , 




(4.2) 


With a c =0, we obtain the simple form, 



| (4.3) 


and Aj = mjj+ — mh r Eq. (4.3) is in agreement with Ref. [68] in the limit a b = 0. 
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Figure 15. 95% Confidence Level allowed regions (white) from fits to the oblique parameters in the 
CP violating 2HDM. 

The result for AS also takes a simple form in the alignment limit [65], 


AS = jcos 2 20 w G(m 2 H+ , m 2 H+ , M%) + sin 2 a b cos 2 a c G(m\ v m\ 2 , M%) 

+ sin 2 a c G(m 2 hl , m 2 3 , M§) + sin 2 a c G(m 2 h2 , M§) + cos 2 a c G(m 2 h3 , M§) 


+ cos ab 


M%) + G{m\ 2 ,m\ 3 , M§) 


+ In 


m hJ n hi rn h 3 


m 


H+ 


Mz) + In [ 


M H,ref 


H+ 


(4.4) 


Analytic results for G(x,y,z ) and G(x,y) are given in the appendix of Ref. [65]. 


- 23 - 

















We use the Gfitter fit to the electroweak data [69], 

S = 0.05 ±0.11 
T = 0.09 ±0.13 

U = 0.01 ±0.11, (4.5) 


with a reference value for the SM Higgs mass, Mfj X ef = 125 GeV. The STU correlation 
matrix is, 


Pij 


( 1 0.90 -0.59\ 

0.09 1 -0.83 , 

y-0.59 -0.83 1 ) 


(4.6) 


and the x 2 is defined as 


A X 2 = Eij(AXi - AX i )(a 2 )~ j 1 (AXj - AXj) , (4.7) 


where Xi = AS, AT, and AU are the central values of the fit in Eq. (4.5), X. L = AS, AT, 
and AU are the parameters in the 2HDM (Eqs. (4.1) and (4.4)), cq are the errors given in 
Eq. (4.5) and o% = cJiPijCFj. 

In Fig. 15 we show the 95% confidence level allowed regions for a & = a c and a c = 0. 
For a/, close to 1, there is some interesting structure due to the interplay of the AS and AT 
limits. For | sina;, |< 0.5, the results are well approximated by the limit from AT only, 


- 80 GeV < A 2 < 120 GeV. 


(4.8) 


4.2 Limits from muon g — 2 

The experimentally measured value of - 9 = a^ places a weak constraint on the parameters 

of the CP violating 2HDM. The deviation between the experimental number and the SM 
theory prediction is [70], 

Aa M = a® xp - a^ M = 265(85) x 10' 11 . (4.9) 

The one-loop contributions from the Higgs sector in the 2HDM to Aare numerically 
small. The larger Higgs sector contributions come from the 2-loop Barr-Zee type diagrams 
with a closed fermion/gauge-boson/heavy-Higgs loop. This class of diagrams can be enhanced 
by factors of M 2 /m 2 ^ relative to the 1-loop diagrams, where M is a heavy Higgs or heavy 
fermion mass. For completeness, these results are given in Appendix C. 

In Figs. 16 and 17, we show the contributions to Aa^ in the 2HDM for relatively heavy 
771 , 2,3 an d rnn+ in units of 1CF 11 . For | sinaj, |< 0.5, there is almost no sensitivity to the CP 
violating phase. The largest contribution is found in the Type-II model for large tan fj and 
is of opposite sign to that needed to explain the discrepancy of Eq. (4.9). 
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Contours of Aa,j in 10 u , type-I 2HDM, a c =ai, Contours of Aa^, in 10 n , type-I 2HDM, ar c =0 




Figure 16. Contributions to (g — 2) p in the CP violating Type-I 2HDM from the Barr-Zee diagrams. 
The heavy scalar masses are fixed to be rrih 2 = 400 GeV, mh 3 = 450 GeV, m H + = 420 GeV, and the 
model parameter v = 1. 
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Figure 17. Similar to Fig. 16 but for Type-II 2HDM. 
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5 Conclusion 


The CP mixture of the 126 GeV Higgs boson is an important property of the Higgs sector 
that deserves further scrutiny. A non-zero CP component is theoretically well motivated and 
may be the origin of the cosmic baryon asymmetry. An important consequence of the 126 
GeV Higgs boson having a sizable CP odd mixture is that the new physics responsible for 
this cannot be decoupled and must lie near the electroweak scale. 

In the context of CP violating, flavor conserving two-Higgs-doublet models, we studied 
the impact of the heavy Higgs searches at the LHC on the CP violating parameters. In this 
class of models, CP violation appears in the neutral Higgs sector, where there are two more 
real scalars (^- 2 , 3 ) in addition to the lightest 126 GeV one. The couplings of the heavy Higgs 
scalars with electroweak gauge bosons are very sensitive to the CP violation in the Higgs 
sector. Turning on a CP odd mixture in the 126 GeV Higgs boson will also turn on the 
heavy Higgs decay channels into gauge bosons, ^ 2,3 —>• WW/ZZ and Zh\. There is data 
from the LHC from the search for a SM like Higgs boson in these decay channels, and the 
non-discovery of a heavy Higgs can be re-interpreted as constraints on the allowed deviation 
from the alignment limit in the two-Higgs-doublet models without CP violation. 

In this work, we point out that heavy Higgs searches are also extremely useful for con¬ 
straining Higgs sector CP violation and in particular the CP mixture of the 126 GeV Higgs 
boson. We demonstrate that the constraints from heavy Higgs searches are largely compli¬ 
mentary to the low energy EDM constraints. We compare our results with the limits from 
the global fit to the 126 GeV Higgs data, and find they can place much stronger limits than 
the light Higgs coupling fit, especially in the interesting regions when there are destructive 
contributions to the EDM. We find in these regions that the heavy Higgs searches are at the 
frontier of probing Higgs sector CP violation. The current limit on the CP violating mixing 
angle, parametrized by 05 , is constrained to be less than 30%, and the LHC heavy Higgs 
search can further narrow down the angle to less than a 10 % level with the high luminosity 
runs. We also expect our work to be a roadmap for the future searches for Higgs sector CP 
violation and the exciting interplay across various experimental frontiers. 

For completeness, we have also explored other relevant constraints from electroweak 
oblique parameters, the muon g — 2 and from B physics, and discussed their implications 
on the heavy Higgs parameter limits in CP violating 2HDMs. 



A Solving the potential parameters in the approximate Z 2 case 


In this section, we list the relations between the potential parameters and the phenomeno¬ 
logical parameters listed in Eq. (2.3) in the approximate Z 2 symmetric 2HDMs. 

rn\ ! = Air ; 2 cos 2 /3 + (A 3 + \±)v 2 sin 2 /3 — Re(m 2 2 e^) tan / 3 + Re(A 5 e 2 ^)r ; 2 sin 2 [3 , (A.l) 
m 22 = A 2 v 2 sin 2 j3 + (A 3 + \a)v 2 cos 2 /3 — Re(m 2 2 e^) cot j3 + Re(A, 5 e 2 ^)t > 2 cos 2 (3 , (A.2) 
Im(m 2 2 ) = r; 2 sin^cos/ 3 Im(A 5 ) , (A.3) 
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B Tri-linear Higgs Couplings 


From the quartic terms in the scalar potential Eq. (2.1), we can obtain the interactions 
between three neutral scalars, in the basis of (H ®, i7 2 ,A°), 


Cis = - (A 0 ) 3 cos/3 ^ 2 sin/3ImAs — cos/ 3 I 111 A 7 


+ -(A 0 ) 2 ■( — 5-ffi cos/3 + Hf cos(3/3) — Ho (5sin/3 + sin(3/3)) 
8 


ReAs 


+4 


i?i cos /3 sin 2 /3Ai + hf 3 cos2 P sin /3A 2 + (Hi cos 3 /3 + iT? sin 3 /3) (A 3 + A 4 ) 

ImA 5 + H% (2H$ - H$ cos(2/3) + H$ sin(2/3)) ImA 7 


+ jA°<j + ((if ?) 2 + (ff 2 0 ) 2 ) sin(2/3) 


+ \ \ H 2 sin/3 


(ff?) 2 A 2 + (Hi) 2 (A 3 + A 4 + ReA 5 ) 


if? cos/3 


(ff?) 2 Ai + (if ?) 2 (A 3 + A 4 + ReA s ) 


(BA) 


- 27 - 





















From these terms one can readily obtain the hihjhk interactions in the mass eigenstate basis 
{h\, h, 2 i / 13 ) using the orthogonal matrix R from Eq. (2.4). In particular, the gm(i = 2,3) 
coefficients used in Eq. (2.24) are 
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d 3 £, 
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(B.2) 


where {H a } = (H ®, H ®, H 0 ). 


C Formula for g — 2 


The magnetic and electric dipole moments of a fermion / correspond to the real and imaginary 
parts of the Wilson coefficient c of the effective operator 


C eff = cfL<J llu f R F ,,v + h.c. , 


(C.l) 


where in the Type-I and Type-II 2HDMs we consider the main contributions to the coefficient 
c that arise from the two-loop Barr-Zee type diagrams. It is straightforward to translate the 
electron EDM results to the corresponding muon anomalous dipole moment. The prescription 
for the translation is, 
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(C.2) 


where AB'y corresponds to those Barr-Zee diagrams with h\ lines connected to the upper 
loop, and the S/F in the bracket corresponds to heavy Higgs scalars/SM fermions running 
in the upper loop. The / 177 , hZ'y and W^H^'y diagram (S) contributions to the EDM have 
been summarized in Refs. [10, 71]. The W^H^'y diagram (F) contributions to the EDM 
vanish in 2HDMs with approximate Z 2 symmetry, but have been calculated in a more general 
framework in Ref. [72]. We perform the above translation based on results in Ref. [72]. See 
also Ref. [73] for a recent work on g — 2 in a 2HDM. 
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We list below the analytic results for the contributions to the muon g — 2 in a 2HDM: 
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The relevant loop functions are, 
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